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Abstract. 

We develop a hydrodynamic description of the classical Calogero-Sutherland liquid: 

a Calogero-Sutherland model with an infinite number of particles and a non- vanishing 
density of particles. The hydrodynamic equations, being written for the density and 
velocity fields of the liquid, are shown to be a bidirectional analogue of Benjamin- 
Ono equation. The latter is known to describe internal waves of deep stratified fiuids. 
Wc show that the bidirectional Benjamin-Ono equation appears as a real reduction of 
the modified KP hierarchy. We derive the Chiral Non-linear Equation which appears 
as a chiral reduction of the bidirectional equation. The conventional Benjamin-Ono 
equation is a degeneration of the Chiral Non-Linear Equation at large density. We 
construct multi-phase solutions of the bidirectional Benjamin-Ono equations and of 
the Chiral Non-Linear equations. 
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1. Introduction 

The Calogero-Sutherland model (CSM) [H [2] describes particles moving on a circle and 
interacting through an inverse sin-square potential. The Hamiltonian of the model reads 



N ^ ^ N 



where Xj are coordinates of particles, pj are their momenta, and g is the coupling 
constant. We took the mass of the particles to be unity. The momenta pj and 
coordinates Xj are canonically conjugate variables. 

The model (classical and quantum) occupies an exceptional place in physics and 
mathematics and has been studied extensively. It is completely integrable. Its solutions 
can be written down explicitly as finite dimensional determinants (for review see [3]). 

In the limit of a large period L ^ oo the CSM degenerates to its rational version - 
Calogero (aka Calogero-Moser) model (CM) where the pair-particle interaction is 
[U The CSM itself is a degeneration of the elliptic Calogero model, where the pair 
particle interaction is given by the Weierstrass p-function of the distance. In this paper 
we discuss the classical trigonometric model ([T]) commenting on the rational limit when 
appropriate. 

We are interested in describing a Calogero-Sutherland liquid^ i.e., the system ([T]) in 
thermodynamic limit when N ^ oo and L ^ oo while the average density N/ L is kept 
constant. We assume that the limit exists and that in this limit a microscopic density 
and current fields 

N 

p{x,t) = ^5(x-x,(t)), (2) 

AT 

jix,^) = ^Pjit)Six-Xjit)) (3) 
i=i 

are smooth single-valued real periodic functions with a period L equal to the period of 
the potential [§] In this case the system will be described by hydrodynamic equations 
written on the density field p{x,t) and the velocity field v{x,t). The velocity is defined 
as j = pv. 

The hydrodynamic approach is a powerful tool to study the evolution of smooth 
features with typical size much larger than the inter-particle distance. Apart from 
application to the CSM, the hydrodynamic equations obtained in this paper are 
interesting integrable equations. We show that they are new real reductions of the 
modified Kadomtzev-Petviashvili equation (MKPl). 

I In the rational case one usually adds a harmonic potential, ^o;^ J2i ^1; to the Hamiltonian to prevent 
particles from escaping. This addition does not destroy the integrability of the system f2^ . 
§ It is likely that there are classes of solutions of the CSM, whose thermodynamic limit consists of a 
number of interacting liquids. In this case the microscopic density give rises to a number of functions 
in the continuum - the densities of the distinct interacting liquids. In this paper we consider a class of 
solutions which leads to a single liquid. 
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In this paper we consider a classical system, however the approach developed below 
can be extended to the quantum case {pj, Xk} = Sjk — > [pj, Xk] = ihSjk almost without 
changes. For a brief description of the hydrodynamics of the quantum system see Ref. |1] . 
The hydrodynamics of the quantum Calogero model has been studied previously O [6] 
in the framework of the collective field theory and some of the results below can be 
obtained in a classical limit (see [7]) of the quantum counterparts of Refs. [51 16]. 

The outline of this paper is the following. In Sec. [2] we parameterize the particles 
of CSM as poles of auxiliary complex fields so that the motion of particles is encoded 
by evolution equations for fields. In Sec. [3] we derive a hydrodynamic limit of these 
equations - continuity and Euler equations with a particular form of specific enthalpy. 
We will refer to these equations as to the bidirectional Benjamin-Ono equation or 2B0. 
We present the Hamiltonian form of 2B0 in Sec. |4} In Sec. [5] we discuss the bilinear 
form of 2B0 and its relation to MKPl. In Sec. M we obtain the Chiral Non-Linear 
equation (CNL) - chiral reduction of 2B0 and discuss some of its properties. In Sec. [T] 
we construct multi-phase and multi-soliton solutions of 2B0 and CNL as a real reduction 
of MKPl. These solutions correspond to collective excitations of the original many-body 
system. Some technical points are relegated to the appendices. 



2. Particles as poles of meromorphic functions 

The Equations of motion of the CSM are readily obtained from the Hamiltonian ([T]) 

Xj=pj, (4) 

V,= -9^ L (^cot-(x,-x.)j . (5) 
^ fc=i (k+j) 

We rewrite this system in an equivalent way as 



■' ^ \ fe=l k=l{kj^j) 




^ 1/)- — I 



N 

.iij g f 'In \ ^ I Uj + Wk \^ Uj + 

Uj 2 \ L J \ ^ — ' Uj — Wk — ' Uj — Uk 



where Wj(t) = e* i""^^*-*^ are complex coordinates lying on a unit circle, while Uj{t) = 
e*i7%(*) are auxiliary coordinates. Indeed, differentiating km with respect to time and 



using (6][7) to remove first derivatives in time one obtains equations equivalent to ( ifS ). 



We note that while the coordinates Xj are real, i.e., \wj\ = 1, the auxiliary 
coordinates, yj{t), are necessarily complex. Given initial data as real positions and 
velocities Xj{0) and Xj{0) one can find complex yj from ^ and then initial complex 
velocities yj{0) from ([7|. Once Xj and xj are chosen to be real they will stay real at 
later times, even though coordinates yi are moving in a complex plane. 
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The coordinates Wj{t) and Uj{t) determine an evolution of two functions 



TV 



9 



TC X - W + Wj 



E 



L '■ — ' w 

TV 



W 



3 



TV 
TV 



TT TT 

— cot — [x 



/I V— \ w -r . \ ^ 

^tZ^ ^ = ^^2^ 



TT TT 

y cot — [x 



Xj), 



w 



w 



e " 



e 



(8) 
(9) 



J i=i 

The latter functions play a major role in our approach. These are rational functions of 
w regular at infinity and having particle coordinates as simple poles with equal residues 
27rg/L. 

The condition that the coordinates of particles Xj are real yields Schwarz reflection 
condition for the function Ui with respect to the unit circle 



Ui{w) = —Ui{l/w) or Ui{x) = — (10) 

where bar denotes complex conjugation. The values of Ui{w) in the interior and exterior 
of a unit circle are related by Schwarz reflection. 

Comparing ([6]), Q and ^ we notice that while the function Ui{w) encodes the 
positions of particles wj, the function uo{w) encodes the momenta of particles as its 
values at particle positions Wj 



TV 

uo{wj)+9j Yl 



L ^ Wj — Wk 

k=i (k^j) ^ 



We notice here that the positions of the particles fully determine the imaginary part of 

(12) 



the field uq on a unit circle. Indeed, we have from (11) 



Im Uo{Xj 



-COt-{Xj 



We now introduce complex functions 

U = Uq + Ml, U = Uq — Ui. 

One can show that they obey the equation 



ut + 



^ 2 , -9 ~ 

-U +t-OxU 



0. 



(13) 



(14) 



Indeed substituting the pole ansatz (8]|9) into (14) and comparing the residues at poles 
Wj and Uj one arrives at ( 6|[7 ). 

The equation (14) connects two complex functions uq and ui. The equation is 
equivalent to the modified Kadomtzev-Petvisashvili equation (or simply MKPl). We 
will discuss its relation to MKPl in Sec. HI 

However, being complemented by the Schwarz reflection condition (10), analyticity 
requirements, and an additional reality requirement it becomes an equation uniquely 
determining uq and Ui through their initial data. 

The analyticity requirements read: Uq{w) is analytic in a neighborhood of a unit 
circle \w\ = 1, while ui is analytic inside \w\ < 1 and outside \w\ > 1 of the unit circle, 
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approaching a constant at w ^ oo. An additional reality requirement is the relation 
between the imaginary part of Uq on a unit circle and Ui stemming from the condition 
(12). We formulate and discuss these conditions in Sec. 3.3 and Sec. |5} 

We will refer to the equation (14) as the bidirectional Benjamin-Ono equation 
(2B0). It is a bidirectional (having both right and left moving waves) generalization of 
the conventional Benjamin-Ono equation (BO) arising in the hydrodynamics of stratified 
fluids [H]. We discuss its hydrodynamic form in the next section. 

The solution of ( |14[ ) given by ^ is the GSM many body system with a finite 
number of particles ([T]). Other solutions describe GSM fluids. They are the central 
issue of this paper. 

To conclude this section we make the following comment. The function ui can be 
expressed solely in terms of the microscopic density of particles ([2]) as 

Ui{w) = -i.gl^^-±^p{C). (15) 

J iTTlQ Q — w 

The integral in this formula goes over the unit circle ^ = exp [i^x) . In the following we 
will denote for brevity p(C) as p{x), when ( lies on a unit circle ( = e^^^ . The density 
itself can be obtained as a difference of limiting values of the fleld Ui at the real x (on 
the unit circle). The discontinuity of ui on the unit circle gives a microscopic density 
(|2| of particles 

Ui{x + iQ) — ui{x — iQ) = —27Tgp{x), lmx = 0, 0<Rex<L. (16) 



3. Hydrodynamics of Calogero-Sutherland liquid 



3.1. Density and velocity 

We assume that in the thermodynamic limit N,L —>■ oo, N/L = const the poles of 
the function ui are distributed along the real axis with a smooth density p(x) and 



consider a complex field Ui{w) given by formula (15). Notice that Ui{w) defined by (15) 
is analytic everywhere outside of the real axis of x (everywhere off the unit circle in 
z-plane) approaching a constant as z — oo. It also satisfies the reality condition (10) 



(the density p{x) is real). In the thermodynamic limit the function Ui is not a rational 
function anymore. It is discontinuous across the real axis with the discontinuity related 



to the density of particles by (16). The value of the field Ui{x) on a real axis (on a unit 
circle in z plane) depends on whether one approaches the real axis from above or below 
(unit circle from the interior z ^^i^ix+io) qj~ from the exterior z — > e*"r(^'~*'^)). More 



explicitly, we have from (15) 

ui{x ± iO) = ngiTP + 



(17) 



The superscript H in the second term of (37) denotes the Hilbert transform and is defined 
as (see Appendix A for definitions and some properties of the Hilbert transform) 

dy „, . , IT 



X] 







^ fiy)cot-{y - x). 
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We also assume that in ^ oo limit the complex field uo{w) remains analytic in 
the vicinity of the real axis in x-plane (i.e., in the vicinity of a unit circle in z-plane). 



The 2B0 (14) does not explicitly depend on the number of particles N . It holds 

oo, N/L = const, however solutions describing a 



also in thermodynamic limit N, L 
liquid are not rational functions any longer. 

We can use 2B0 to define velocity through the continuity equation 

Pt + d^ipv) = 0. 



(19) 



The discontinuity of the complex field u{x) (13) across the real cLXlS, clS well as a 



discontinuity of the field ui (see (16)) is the density 

u{x + iO) — u{x — iO) = ui{x + iO) — ui{x — iO) 



-2ngp{x). 



(20) 



Differentiating (20) with respect to time and using 2B0 (14) we obtain the continuity 



equation and identify the velocity field v{x) as 

1 



vix) 



^io(a;) + - {ui{x + iO) +Ui{x-iO)) - igd^hg y/p{x) 



+ W {-^P^ix) - log ^/p{x)) 



or 



Uq[X\ 



- W (ttp - log y/p) 



(21) 



(22) 



Since v{x) is a real field (22) provides a reality condition analogous to (12). Indeed, one 



can see from (22) that 



lmuo{x) = -g (vrp^ - 9^ log ^p) , (23) 

i.e., the imaginary part of Uq{x) is completely determined by the density of particles or 
equivalently by the field ui. It is also convenient to have an expression for u{x) on a 
real axis 



u{x ±iO) = v + g (=F7rp + id^ log ^p) . 
It has the same discontinuity across the real axis as ui{x). 



(24) 



3.2. Hydrodynamic form of 2B0. 



Now we are ready to cast the equation ( 14 ) into hydrodynamic form. 



Taking the real part of 2B0 (14) on the real axis and using identifications ( 17p2 ) 



and the continuity equation, (19), after some algebra we arrive at the Euler equation 

vt + d^,(^^+w{p)^ =0, (25) 
with specific (per particle) enthalpy or chemical potential] given by 



wi.p) = 7;i^9p) 



1. g' I 



dly/p + ng^p^ 



2^H 
X ■ 



2^-"' 2Vp 

The specific enthalpy and chemical potential are identical at zero temperature. 



(26) 
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Equations ( I9p5 ) are the continuity and Eulei{^ equations of classical Calogero- 
Sutherland model. They are the classical analogues of quantum hydrodynamic equations 
that have been obtained for the quantum CSM in Refs. |5l El |9] first using collective 
field theory approach [TOl [TTl [12] and later by the pole ansatz similar to the one used 
above It was noticed in [T2j and then in [7] that the system ( p!9|25[26[ ) has a lot of 
similarities with classical Benjamin-Ono equation [13]. The similarities and differences 
with Benjamin-Ono equation are discussed below. We will refer to ( l9|25[26[ ) as to a 
hydrodynamic form of the bidirectional Benjamin-Ono equation (2B0). 



3.3. Bidirectional Benjamin-Ono equation (2B0). 
Let us now summarize the 2B0 equation: 



^ 2 , -^o ~ 

-U + 1-On.U 

2 2 



0, 



Ut + 9^ 
U = Uq + Ml, M = Mo — Ml. 

The functions mq and mi are subject to analyticity conditions 

Ml (a;) - analytic for Im (a;) ^ 0, 

Uq{x) - analytic for |Im(x)| < e for some e > 0, 

and to reality conditions 

Mi(a;) = —ui{x). 



(27) 
(28) 

(29) 
(30) 

(31) 



In addition, the fact that the equation (27) holds in the upper half plane and in the 



lower half plane (inside and outside of the unit circle) yields the condition 
Im [u{x ± zO)] = log Re [Mi(a; ± zO)]. 



(32) 



It also follows from ( [I7p3p4l ). The condition (|32|) looks more "natural" in the bilinear 
formulation (see eq. ([56]) below). 

These reality and analyticity conditions reduce two complex fields mq and ui to two 
real fields - density p{x) and velocity v{x) as ( 17p2 ). Then, a complex equation ( [l4| ) 
defined in both half planes immediately yields the hydrodynamic equations ( [T9|25[|26| . 
Inversely, knowing real periodic fields p(x) and v{x) one can find fields mq, mi everywhere 
in a complex x-plane. 



Mode expansion The analyticity and reality conditions can be recast in the language 



of mode expansions. It follows from (15) that 



-ng (po + 2 J2n=i PuW"") 
irg (po + 2 E^^i p^w-") 



\w\ < 1 
\w\ > 1 



(33) 



^ The eq. ( 25 1 has a form of an Euler equation for an isentropic flow. Because of the long range character 
of interactions the enthalpy cannot be replaced by the conventional pressure term dxw{p) — > p~^dx{pip)) 
- the standard form of the Euler equation. 
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where p„ = pl„ = jj^ ^ p(x)e are Fourier components of the density. 

The values of the field Ui{w) in the upper and lower half-planes are then 



automatically related by Schwarz reflection (10). 



Conversely, the field uo{x) being analytic in a strip around the unit circle is 
represented by Laurent series 



K) + ^ (flnU^" + hnW ") , |Im logw| < 2'Ke/L. 



(34) 



n=l 



The 2B0 equation remains intact in the case of rational degeneration. Rational 
degeneration of formulas of the Sec. [2] are obtained by a direct expansion in 1/L . In this 
limit fields are defined microscopically as Mi(x) = —ig Y^- ^^73^ and Uq{x) = ig J2j ^r^- 



4. Hamiltonian form of 2BO 



The 2B0 is a Hamiltonian equation. Let us start with its Hamiltonian formulation in 
the hydrodynamic form pt = {H, p} , Vt = {H, v} with the canonical Poisson bracket of 
density and velocity fields 

{p{x),v{y)} = 5'{x-y). (35) 

Equations ( p!9|25[26[ ) follow from 



H 



9^ 



dx { ^ + pe{p) 



(vrp^ - 9^ log ^f. 



(36) 



(37) 



Here the "internal energy" (37) and the enthalpy (26) are related by a general formula 

For references we will give alternative expressions for the Hamiltonian. Let 
\E' = ^/pe^^ where v = gdx^ then 

-.2 



H 



9 



\dx^ - T[p^'^\ dx, 



(38) 



where p = l^&p. The Poisson's brackets for \E'(x) are canonical: {\E'(x), \E'(y)} = 0, and 



{\E'(x), \E'*(y)} = — y). The equations of motion for ^ and are 



9 



(39) 



and its complex conjugate. A simple change of a dependent variable $ = vE^e*'^-^ l**-^' 
leads to 



9 



(40) 



where /"'" denotes the function analytical in the upper half-plane of x defined as 
/+ = One can recognize in MOh the intermediate nonlinear Schrodinger equation 



(INLS) which appeared in Ref.fT^ as an evolution of the modulated internal wave in a 
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deep stratified fluid. Therefore, one can alternatively tliink of 2B0 as the hydrodynamic 



form of ( |40[ ) identifying hydrodynamic flelds p and v to be 

i 



$ = i/pexp 



9 



dx {v + Tigp) 



or with the fleld u{x) from (24) as 
igd:^ log $ 



MX 



iO). 



The Hamiltonian (|36|) or (|38|) can be rewritten in terms of $ as 
H 



2 



9^.$ - ■i27rp+$ dx, 



(41) 



(42) 



(43) 



where p = However, the Poisson's brackets for $ are not canonical anymore]^ 

2B0 is an integrable system. It has inflnitely many integrals of motion. The 
flrst three of them follow from global symmetries. They are conventional the number 
of particles = J dx p, the total momentum P = J dx pv, and the total energy 
H = J dx + pe{p)j. They are conveniently written in terms of the flelds u and 
u as 

1 



h 



N 



dxu, 



dx -u^, 



2vr^ Jc 

-/ 

27r^ Jc 

= (p dx 

27r^ Jc 

where the integral is taken over the both sides of the unit circle, ("double" contour C 
shown in Fig. [T]). For more details on conserved integrals see Appendix B 



2H 



^ 3 , o ~ 

-U + l-UOxU 

3 2 



(44) 
(45) 
(46) 



The Poisson bracket for the flelds Mo(w) and Ui{w) can be easily obtained from 
from ( 17p2 ) and (35) by analytic continuation. We flnd that {uq{w),Uo{w')} = 
{ui{w) , ui{w')} = and 



{uoiw),Ui{w')} = ig 



WW 



L J {w — w') 



i\2 



igd^^ cot ^{x-y). 



(47) 



5. Bilinearization and relation to MKPl equation 



The equations described in the previous section, their integrable structures and their 
connection to integrable hierarchies are the most transparent in bilinear form. 
Let us introduce tau-functions tq and ri as 



Mo 

Ml 



igdx^ogTo, 
-igdx log Ti. 



(48) 



Simple calculation using (35l gives {^{x),^{y)} = ^<i>(2;)<i>(?/) sgn (x - y), {(f>(a;), <&*(?/)} = lS{x - 



y) — ^<^(x)<^*{ii) sgn [x — y) and similar expressions for complex conjugated fields. One should tliink of 
'^{x) as a canonical bosonic field while of ^{x) as a classical analogue of a field with fractional statistics. 
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Figure 1. Contour C surrounding the unit circle is shown together with our 
conventions in defining right and left fields. 



It can be easily checked that the 2B0 (14) can be rewritten as an elegant bilinear Hirota 
equation on r-functions: 



0. 



Here we used the Hirota derivative symbols defined as 
Dlf{x) ■ g{x) 



For example, 



Dtf-9 

Dlf-g 



\iYii{d^-dyYf{x)g{y). 



idtf)g-fidtg), 
{dlf)g-2{dj){d,g) + f{dlg), 



(49) 



(50) 



(51) 



etc. 



We emphasize that the bilinear equation holds on both sides of the unit circle. 
Introducing notations 

r±i=ri(x±iO) (52) 

we can rewrite the equation as 



2 ^ 



T-1 ■ ro= 



0, 



9 



-Dl 
2 



^0 ■ r_i 



0. 



(53) 
(54) 



Equation (49) is the modified Kadomtsev-Petviashvili equation (MKPl). MKPl 
contains two independent functions ti and tq and is formally not closed. The analyticity 
and reality conditions (29 32| stemming from the fact that all solutions are determined 
by two real functions p{x,t) and v{x,t), close the equation. Under these conditions the 
equations can be seen as a real reduction of MKPl. Let us formulate these conditions 
in terms of tau-functions. 
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The first requirement is that t±i is analytic and does not have zeros for Ima; > 
(< 0) after analytic continuation. Also tq should be analytic and should not have zeros 
in the vicinity of the real axis, i.e., for |Imx| < e for some e > 0. 



The second requirement is that t±i should be related by Schwarz reflection (10). 
In terms of tau-functions it becomes on the unit circle (for real x) 

r_i=r^Te^®W, (55) 

where a phase 0(t) can be any time-dependent function. 

The third requirement is related to the fact that Im uq is a function of density only 
and, therefore, can be expressed in terms of ui as can be easily seen from ( [I7p2| ). This 



condition ( 32 ) can be written in a bilinear form as follows 



(56) 



The multiplicative constant in the r.h.s of (56) fixes the relative normalization of tq and 



Ti and is arbitrary. We have chosen it to be 27r. 

Finally, we note that the pole ansatz solution (8]|9) corresponds to the polynomial 
form of tau-functions with zeros at Wj and Uj 

N 
N 



(57) 



(5J 



i=i 



6. Chiral Fields and Chiral Reduction 



6.1. Chiral fields and currents 

The 2B0 equation can be conveniently expressed through yet another right and left 
handed chiral fields 



Jr^L = V±g[Kp + 9^. (log y/p)^] . 

These fields are real. In terms of them, the 2B0 equation (14) reads 



(59) 



dtJR.L + 



Tgd. [jR,Ld.{\og^p)'' - (J«,L9,logv/p)^] = 0. 



(60) 



Here p is a function of Jr and Jl implicitly given by (59). The Hamiltonian acquires a 
Sugawara-like form 

1 



H 



dxp {jR + JLf + {J^-Jt 



H\2 



(61) 



* Jr.l can be expressed solely in terms of uq field. It is easy to check that (59 1 is equivalent to 
Jr,l = Re (tto T iuo)- 
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with Poisson brackets 



{jR,L{^),jRAy)}= ±27igdJ{x - y) 



± —rdxdy 
2L ^ 



1 



1 



{JBix),JL{y)} 



_9_ 
2L 



dxdy 



1 



pi.y) 

1 



TT 



cot -{x-y) 



p{x) p{y) 



TT 



cot -{x-y) 



(62) 



(63) 



We note that Poisson brackets become canonical and left and right fields decouple in 
the limit of a constant density. 



6.2. Chiral Reduction 



We first note that the right and left currents Jrl are not separated in eq. (60). The 



equations for J/j and are coupled through the density p which should be found in 



terms of Jr^l from (59). However, it is possible to find the chiral reductions of 2B0 



assuming that one of the currents is constant. We explain this reduction in some detail 
in this section. 



The 2B0 (14) or (60) admits an additional reduction to a chiral sector where 
one of the chiral currents (59), say left current, is a constant JL{x,t) = vq — rrgpo. We 
can always choose a coordinate system moving with velocity vo- This is equivalent to 
setting the zero mode of velocity to zero vq = 0. The condition = —irgpo becomes 

v = g [nip - po) + dxihg y^f] . (64) 

Then the currents can be expressed in terms of the density field only 

J^(x) = Jo, 

jRix) = Jo + Jix), (65) 

Jo = T^gpo, 

J{x) =2g[n{p-p^) + dxi\og^pf]. (66) 
It follows from Eq. (60) that once the current Jl is chosen to be constant Jl(x) = Jq at 



t = it remains constant at any later time. The condition (64), therefore, is compatible 



with 2B0. Then the density p{x,t) evolves according to the continuity equation (19) 



with velocity determined by the density according to (64). We obtain an important 



equation (written in the coordinate system moving with velocity vo) 

pt + g[p (vr (p - Po) + dxilog ^f)]^ = 0. (67) 
We refer to this equation as the Non-Linear Chiral Equation (NLC). A substitution of 



the chiral constraint (65) to (61) gives the Hamiltonian for NLC 



H=lldxp[j' + {j"f] 



(68) 



with Poisson brackets for J{x) following from (62). This equation constitutes one of 
major results of this paper. 
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NLC can be written in several useful forms. One of them is: 



where p{x) = po e' 



ip{x) 



vrpo (2e'^ + -y^f 



, 9 H 



(69) 



6. 3. Holomorphic Chiral field 



Under the chiral condition (|64|) the field uq becomes analytic inside the disk. Indeed 



combining (64) and (22) we obtain 

d( C + w 



Uo[W 



2TTi( ( — w 



\w\ < 1. 



(70) 



In the chiral case it has only non- negative powers of w in the expansion (34). Negative 



modes vanish 6„ = 0. Conversely, the condition of uq to be analytic inside the unit disk 
is equivalent to Jl = const. 



The current itself (65) is the boundary value of the field Rewo harmonic inside the 



disk 



2 Jo + 2ReMo = 2Jo + {anw'' + a„w"" 



) 



(71) 



n=l= 



The fields u and u are in turn also analytic inside the disk. Let <y9 to be a harmonic 
function inside the disk with the boundary value log(p/po)- Then (p = (f){w) + (t){w), 
where (j){w) = {\og{p/ po))~^ . Here f~^{w) = | / ^r^/(C) 2^ ^ function analytic in the 
interior of a unit circle which value on the boundary of the disk is (/(x) — if^{x))/2. 
It follows from pipi) ) that 

-Jo + igdcj), \w\ < 1, (72) 

M + 47r(7po(e^)+, \w\ < 1, (73) 



u 



u 



Then 2B0 (27) becomes an equation on an analytic function in the interior of a unit 
circle 



+ i^-[{d<pf + d^(f\+T^gpod{2e 



0. 



(74) 



This is the "positive part" of (69) which is a direct consequence of (67). 



We remark here that the chiral equation (67) has a geometric interpretation as 



an evolution equation describing the dynamics of a contour on a plane. Within this 
interpretation the term ^^^(log ^/p)^ of (67) is the curvature of the contour (see Appendix 
C 



6.4- Benjamin- Ono Equation 



Another form of the Chiral Equation (67) arises when one considers the fields u and 



u outside the disk. There neither u nor u are analytic, but their boundary values are 
connected by the Hilbert transform 

u{x -tO)= - Jo + 2g [np + ^^^.(log v^)+] , (75) 
u{x -iO) = - Jo- iu"{x - iO). (76) 
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The bidirectional equation Eq 
unidirectional (chiral) 

~1 



(27) complemented by this condition becomes 



Ut + 5a; 



2 , 9 ^ 

-u H — OxU 
2 2 



H 



0. 



(77) 



This is just another form of the chiral equation (67) 



The chiral equation (77) has the form of the Benjamin-Ono equation [13]. There 



are noticeable differences, however. Contrary to the Benjamin-Ono equation, Eq. (77) 



is written on a complex function, whose real and imaginary values at real x are related 



by conditions ( 75 ) implementing the reality of the density: 

'Reu = -jQ + 2gTxp + g{dx\og^/p)^ , Imu = gd^log ^/p. (78) 

One understands this relation as a condition on the initial data. Once it is imposed 
by choosing the initial data for the density p, the condition remains intact during the 
evolution. 

However, in the case when the deviation of a density is small with respect to the 
average density \p — po\ po, the imaginary part of u vanishes in the leading order of 
1/po expansion 

Jo + 2ngLp ^ 2ng{p - po) + -''o 



and the condition (78) becomes non-restrictive. In this limit Eq. (77) becomes an 



equation on a single real function. It is the conventional Benjamin-Ono equation. 



One can think of NLC (67) as of finite amplitude extension of BO. Similarly, 2B0 



is an integrable bidirectional finite amplitude extension of BO. It is interesting that 
there exists another bidirectional finite amplitude extension of BO - the Choi-Camassa 
equation ^6\. However, it seems that the latter is not integrable. 



7. Multi-phase solution 

In this section we describe the most general finite dimensional solutions of 2B0. These 
are multi-phase solutions and their degenerations - multi-soliton solutions. In the former 
case the r-functions are polynomials of e*'^*'^', where ki is a finite set of parameters, 
the latter are just polynomials of x. These solutions are given by determinants of 
finite dimensional matrices. They appeared in the arXiv version of Ref . [T7] . One can 



construct those solutions using the transformation (41) of 2B0 to INLS (40). For the 



latter multi-phase solutions were written in [H] (see also [IHl [19]). We use a different 
route in this section deriving multi-phase and multi-soliton solutions as a real reduction 
of corresponding solutions for MKPl. 

7.1. Multi-phase and multi-soliton solutions of MKPl 

We start from a general multi-phase solution of MKPl equation and then restrict it to 
2B0 equation. 
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A general multi-phase solution of MKPl equation 
'zA + ?/^'Vi-ro = 0. 



is given by the following determinant formulae 



det 



^jk ~l~ ^a,j 



Pj - Qk 



0,1 



Cij _ qj- K - Vq 

Co J Pj- K -Vq 

where the phases are 



gOj{x,t) = {qj -Pj){x - Xoj) 



Pi 



t. 



g 9o{x, t) = Kx-^t- Uvo + K)x - ^^^^±^t 



g9i{x,t) 



Kx 1. 

2 



(79) 

(80) 
(81) 

(82) 
(83) 

(84) 



This solution is characterized by an integer number (number of "phases"), and by 
AN — 1 parameters p-,-, g-,-, cqj, x^j and moduli K and vq. The solutions become single- 
valued on a unit circle if pj and q^ are integers in units of g^. 

7.2. Multi-phase solution of 2B0 



Without further restrictions the parameters entering (80 84) are general complex 



numbers. Reality nature of 2B0 equation restricts them to be real. 

The real moduli K and vo are obviously zero modes of the fields ui and uo 
respectively, and therefore, they are zero modes of the density po = f pdx = —K/ (ng) 
and velocity J v dx = vq. 



7.2.1. Schwarz reflection condition We have to restrict the coefficients Caj, so that 



there exists another solution r_i,ro of Eq. (79) sharing the same Tq with the solution 



(86|87) and obeying the Schwarz reflection property (55). 



The Galilean symmetry of the equation (79) is here to help. If Ta{x,t), a = 0, 1 
give a solution of (79) then the pair e*'^''^~*^'"*ra(x — gPat,t) is also a solution provided 



that Ea+l-Ea = UPa+l-Pa 



Being applied to the solution (80 84) the Galilean invariance can be utilized as 



follows. We notice from (82) that 



9j{x - Pt, t- fe, qj}) = Ojix, t- {pj + P, qj + P}). 



^5) 



Performing the Galilean boost to (80), multiplying both tau-functions by e and 
shifting Pj — > Pj — P, qj qj — P we obtain that 



jei~iP{K+vo)t+^{Px- i^t) 



det 



5jk + 



Pj -qk 



qj - K -Vq- P 
Pj - K -Vq- P 



(86) 
(87) 
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form a solution of (79) with the same tq (80) 



Now we are going to show that for a particular choice of coefficients Ci j the Galilean 
boosted solution ( 86p7 ) is a complex conjugate of r+i from (80). To show this we will 
employ the determinant identity (D.7)j|] 

We apply the determinant identity ( |D.7[ ) to ( [86| and obtain 



T-l 



jei ~iP{K+vo)t+ ^ (Px-^^t) 



n 



det 



6jk + bj 



Pj -qkj 



where 



The Schwarz reflection condition(|55|) r_i = r+ie*®*^*^ requires 
bj = cij, 

determines Q{t), and gives a relation 



(89) 

(90) 
(91) 



Finally, combining all relations ( [8l|87[89|[90| ) together we obtain the condition on 
coefficients Caj 



Pi - Qj 



n 



{Pj - qk){qj - Pk) 



(92) 



Qj- K - Vo 



Pj- K - ^;o V "^^ Pj - K - Vo - P 



Qj - K -Vq- P' 



Condition (92) is necessary to turn a general solution of MKPl into a solution of 



2B0. We also have to find a condition that guarantees that ti has no zeros inside the 



unit disk. Before turning to this analysis, we first discuss degeneration of formulas (80) 
into multi-soliton solution. 



7.2.2. Multi-soliton solution of 2B0 The multi-soliton solution of 2B0 follows from 
the multi-phase solution in the limit pj qj. We introduce 

kj = Pj - qj, (93) 



= T^iPj + Qj) 



(94) 



' 2' 

and consider the limit kj keeping Vj fixed. After some straightforward calculations 
we obtain 



Ta 



z'^^ det 



5jk (x - XQj - Vjt + iAaj) + ig 



6 



'jk 



A, 



9 



a,] 



2 \Vj — vq + K 



± 



1 



'3 

a = 0,1. 



(95) 
(96) 



jj A similar trick was used by Matsuno [TH] to prove the reality of a multi-phase solution for conventional 
Benjamin- Ono equation. 
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One notices that in the hmit t — > -\-oo the solution (95) asymptotically goes to the 
factorized form 

3 

describing separated single solitons. 

Eq. (97) gives a large time value of zeros of ri. Their imaginary part is 

-ReA^^,=g- \ (98) 

It must be negative in order for ri to have no zeroes inside the unit disk. Since A' < 
we must require 

{v,-v,f>K\ (99) 

In the next paragraph we argue that under this condition and additional restrictions on 
parameters pj, qj (see eq. (109) below) the moving zeros never cross the real axis, and 
therefore zeros stay outside of the unit disk at all times. 

To conclude this section we note a unique property of the 2B0 equation (shared 
with the BO equation). Namely, there is a "quantization" of the mass of solitons: each 
soliton of 2B0 carries a unit of mass regardless of its velocity. We have for iV-soliton 
solution 

j dx{p- po) = A^. (100) 

Where K = —ngpo. The total momentum, and the total energy of a multisoliton 
solution is given by 

/ dx {pv - povo) = ^ Vj, (101) 
/ + P<P) - ^ - Poe(po)) = E f ' (102) 



where e(p) is defined in (37). 



One-soliton solution has a form 

1 A, 

IT {x - xoi - vity + Af 
Ao 

{x-xoi -vity + Al 



P = Po + Z7Z 1 ^\'?. I (103) 

1 



where 



A, = Re = -, (105) 



{vj -voY - (irgpo) 



Aq = ReAo,i = 7 ^ r^. (106) 



[vj -VoY - {irgp^ 



'0) 



This one-soliton solution has been found first in Ref. [7] (see also [23] )• 
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7.2.3. Analyticity condition Now we can turn to the multiphase solution and derive 
conditions sufficient in order for Ui to be analytic in the upper half-plane in complex 
X- variable (inside the unit disk). Analyticity in the lower half-plane follows from the 



Schwarz reflection condition (|55|). We will follow the approach of Dobrokhotov and 
Krichever 



developed for Benjamin-Ono equation. 



Analyticity of ui means that ri given by (80) has no zeros in the upper half plane, 
or that the matrix 



M. 



jk 



is non-degenerate. 



- Ojk H 

Pj - Qk 

Following the approach of Ref. 



(107) 



23] we derived in Appendix E 



sufficient condition of non-degeneracy of the matrix M from (107). Let us now write 



that condition ( E.10[|K9| with Cj defined by ( E.l|[92 ). We obtain (calculating fj] 
P T-r Pj - Qk 



n 



where we used shifted numbers pj 

-2K + z,iP,-q,)- 

The set of conditions 



Pk 



same sign for all j, 



108) 



Pj — K — vo and similar for qj. Here P 



qi < pi < . . . < qm < Pm < < P < qm+i < Pm+l < ■ ■ ■ < qN < Pn (109) 



satisfies (108). Moreover, (109) yields to (99), which in its turn means that at least at 



some values of parameters (large time and soliton limit) no zeros of ti are inside the 
unit disk. Since they also can not be on the circle they do not cross it while moving in 
time and in the space of parameters. 



Condition (109) suggests that a general solution is characterized by a integer 



number — 2m. This is chirality - the difference between the number oi N — m 
right and m - left moving modes 
1 



27ig 
1 

2^ 



(Jr -Vo- ngpo) dx = N - m, 



{Jl -Vo + ngpo) dx 



m. 



(110) 



1111 



Eqs. (|80 p4p2|l09| summarize a general finite dimensional quasi-periodic solution. We 
emphasize here that this solution is not chiral and contains both right and left-moving 
modes. 



7.3. Multi-phase solution of the Chiral Non-linear Equation 

The (right) chiral case appears when tq has no zeros outside the unit disk. It naturally 
happens when the number of, say, left-moving modes m in (109) vanishes m = 0. In 
this case all Vj —Vq<0. In their turn, imaginary parts of zeros of tq in the multi-soliton 

>0 (112) 



limit (as in (98)) 
-A 



Oi = -9 



(v, - vor - 
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are positive. One can check that in this case (107) with ci ,,• cqj is non-degenerate for 
arbitrary values of parameters satisfying (109) with m = (and similarly for m = N). 



Therefore, to{x) is non-zero in one of half-planes. 
This is a chiral multi-phase solution of 2B0. 



7.4- Multi-phase solution of the Benjamin- On o equation 

The known solutions of the Benjamin-Ono equation [201 El] are obtained from the 
solutions of the Chiral Non-linear equation by taking the limit po —>■ oo. In this case 
K — oo and conditions (109) allow for a good limit only if m = (left sector) or 
m = (right sector). Let us concentrate on the right sector. We redefine pj pj — K, 
1j ~^ 1j ~ "^0 "^0 ~ 2K, go to the frame moving with velocity —K (x — > x -|- Kt), 
and obtain from ( 109|[92 ) in the limit K ^ —oo 

(113) 
(114) 



qi < Pi < q2 < . . . < Qn < Pn, 

2 



C, 



n 



{Pj -Pk){qj - Qk) ^ / Pj - Vq 
(Pj - Qk)iqj - Pk) \qj-Vo 



l-2a 



with solution given by ( [80p2| ) and with ( [83p4l ) (one should put = in latter two). 
This is nothing else but the multiphase solution of conventional Benjamin-Ono equation 

mm. 



7.5. Moving Poles 



The 2B0 equation ( 14 ) looks very similar to the classical BO equation. One of important 
tools in studying the classical BO equation is the so-called pole ansatz - solutions in 
the form of poles moving in a complex plane. [20] We have already seen that the pole 



ansatz (8]|9) describes the dynamics of the original Calogero-Sutherland model with 
finite number of particles N. 

In this section we consider collective excitations of Calogero-Sutherland model in 
the limit of infinitely many particles. These excitations are given by "complex" pole 
solutions of the 2B0. 

In the Pole Ansatz (8][9), the reality conditions were satisfied by requiring Xj to be 
real (or Wj{t) moving on a unit circle). One could generalize the Pole Ansatz (8][9) to 
case where Wj{t) are away from the unit circle and moving in a complex plane. The 
equations (6][7) describing the motion of poles preserve their form. However, u^i{w) 
outside of the unit circle is not related to the ui{w) inside of the circle by analytic 



continuation but only by Schwarz reflection (10). The field Ui{w) is analytic inside the 



unit circle and has moving poles outside of the unit circle (and vice versa for U-i{w)). 
Of course, having obtained the solution of 2B0 inside the unit circle does not mean 



automatically that the Schwarz reflected function ( 10 ) will solve 2B0 in the exterior of 



the circle with the same uq. The property (10) requires that (6][7) are satisfied not only 
by Uj and Wj but also by Uj and 1/wj. This requirement will significantly constrain the 
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positions of poles Wj and Uj in a complex plane. It turns out that this constraint allows 
for non-trivial solutions. 

We emphasize here once again that while real axis poles Xj oi ui in the pole ansatz 
represent the original CS particles, the complex poles Xj represent collective excitations 
of the CS liquid moving in the background of macroscopic number of particles. 

Instead of looking for moving pole solution in this section we have taken a different 



route. We first construct the much more general solution of 2B0 (14) with proper 



reality conditions and then obtain a moving pole (i.e., multi-soliton) solution as a limit 



of the multi-phase solution. One can see from (95) that for soliton solutions the zeros 



of tau-functions move in a complex plane. It is especially clear at large times when 



solitons are well separated (97) 



8. Conclusion and discussion 

In this paper we have shown that the dynamics of the classical Calogero- Sutherland 
model in the limit of infinite number of particle is equivalent to the bidirectional 



Benjamin-Ono equation (14). The bidirectional Benjamin-Ono equation (2B0) is an 
integrable classical integro-differential equation. Its integrability can be deduced from 
the fact that it is a Hamiltonian reduction of MKPl as it is shown in this paper. As an 



alternative, one can use the equivalence of 2B0 to INLS (40). The integrability of INLS 
was proven and the spectral transform was constructed for INLS in Ref. |25] (see also 
[19]). Therefore, one can use all techniques developed in the field of classical integrable 
equations for 2B0. It has multi-phase solutions (explicitly constructed in this paper), 
bi-Hamiltonian structure, an associated hierarchy of higher order equations, etc. 2B0 
is intrinsically simpler than many other classical integrable models. Its solitons have 
"quantized" area independent of soliton's velocity. The collision of two solitons goes 
without any time delay etc. This is a reflection of the fact that underlying Calogero- 
Sutherland model is essentially a model of non-interacting particles in disguise. In 
particular, 2B0 supports a phenomenon of dispersive shock waves. Some applications 
of this phenomenon to quasi-classical description of quantum systems were considered 
in [ig. 

Most of the results of this paper can be generalized along two avenues: 
generalization to an elliptic case and generalization to a quantum model. 

The Calogero-Sutherland model (trigonometric case) can be generalized to an 
elliptic case — elliptic Calogero model where the interaction between particles is either 
Weierstrass p(x|ti;i, ti;2)-function with purely real and purely imaginary periods uji,iuj2, 
and to its hyperbolic degeneration (hyperbolic case) with inter-particle interaction given 
by sinh~^(a;/a;2) (see [3] for review). 

In both cases most of formulas remain unchanged if one substitutes the Hilbert 
transform for a transform with respect to a strip < Imx < u, where a; is an 
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imaginary period 

= / ((x - or j — coi\i—{x - x')f{x')dx' . (115) 

J J UJ2 UJ2 

In the first case the integration goes over a real period of the Weierstrass (^-function. 

The elliptic Calogero model allows one to study a crossover between liquids with 
long range inter-particle interaction to liquids with short range interaction. In the limit 
of a large imaginary period u;2 — > cxd the p-function degenerates to 1/ sin(a;/u;i)^ - the 
case of long range inter-particle interaction. The opposite limit a;2 — > gives rise to a 
short range interaction: uj2p{x) 5{x). 

In the latter case the the Hilbert transform (115) becomes a derivative ^d^f 
and the equations discussed in this paper become local. In particular the Benjamin-Ono 
equation flows to the KdV equation, while the bidirectional BO-equation flows to NLS 
- the nonlinear Schrodinger equation. 

2B0 in the limit of small amplitudes and in the chiral sector becomes the 
conventional Benjamin-Ono equation. In elliptic case (and in the hyperbolic one) the 
limit of small amplitudes in the chiral sector leads to a generalization of the Benjamin- 
Ono equation, known as the ILW (Intermediate Long Wave) equation [S]. Contrary to 
the Benjamin-Ono equation and to 2B0, the latter and its bidirectional generalization 
2ILW have elliptic solutions. 

We intend to address the elliptic case in a separate publication. 

Probably, even more interesting is a generalization of the results of this paper to 
the quantum case. It is well known that the classical GSM model ([T| can be lifted to a 
quantum integrable Calogero-Sutherland model [HEIET]. The latter model is defined by 
([T]) with — i> h^\{X — 1) and Pi = —ihdx^- The 2B0 equation in the form (14) remains 
unchanged, except for the change of the coefficient g X — 1 and for the change of 
Poisson brackets (47) by a commutator: {, } — > |[, ]• The change g ^ \ — 1 valid 
for eq. (14) is not correct for all formulas. For example, the bilinear form of classical 
2B0 (49) is identical to its quantum version with just a change of notations g X. 
For some details see m. Multi-soliton solution of 2B0 presented here corresponds to 
exact quasiparticle excitations of quantum Calogero-Sutherland model [29l [7] . The more 
detailed study of the relations between integrable structures of the classical 2B0 and 
its quantum analogue is necessary. 
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Appendix A. Hilbert transforms 

Given a function /(x), f{x) as x ^ ±00, tiie Hilbert transform is defined as 

f"{x) = H^^ dy^. (A.l) 

T^J^oo y-x 

For periodic functions with period L we define tlie transform as 

f''{x) = f^fiy)cot^iy-x). (A.2) 



L L 

The Hilbert transform of the constant function is zero. The Hilbert transform is inverse 
to itself = —1 or (f^)^ = —f and it commutes with derivative {f^)x = ifx)^- 
More generally, a function /(x) defined on the closed (directed) contour C surrounded 
the origin of a complex plane can be decomposed into a sum / = + f~ , analytic 
functions inside (outside) of the contour, such that /^(O) = /^(oo). Then 



Using (A. 3) it is easy to derive the following properties: 

f'g + ir = ifgf - (/^^^)^ (A.4) 

and some integration formulas 

dxf^ =0, (A.5) 

dxf''g= - J dx/g"", (A.6) 

dxff = Q. (A7) 

From (A. 3) we have for functions analytic in one of half planes 

{f^f = ±if^ (A.8) 
and as an immediate consequence 

l^^ikx^H =ie*'=^sgnfc, (A.9) 

^ for Ima >0. (A. 10) 



X — a J X — a 

Generally, in Fourier space the Hilbert transform is equivalent to a multiplication by 
isgn/c, i.e., for Fourier coefficients 

(/^)fc = z(sgnA;)/fc. (All) 

It is easy to derive the following useful identities 

u^dx= j{u"fdx, (A.12) 

u^dx =3 / u{u^fdx, (A.13) 
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u dx 



u 



j [{u''f + Qu\u''f] dx, 
dx = h j [u{u''f - 2u\u''f] dx. 



(A. 14) 
(A.15) 



Appendix B. Conserved integrals of 2BO 



In this section we present tlie conserved integrals of 2B0 written in different forms. The 
contour integrals below are taken along the contour C defined in Figure [T] 

C dx f M^; Tl 

(t u = dx p = / — tOx log — . 

J 2ng J ^ y 27r " ^ r„i 

[ dx 1 2 ([ 



h 
h 



(B.l) 



J 2Tig 2 J 2'ng 



UiUo 



g r dx 
2 J 2^ 

/dx 
2^g 



DIt-i ■ To Din ■ To 



T-lTo 

^ 3 , o ~ 

-U + l-UOxU 

3 2 



riTo 



dx pv 



idt 



dx , Tl 

TT — ■ 
2tt r_i 



(B.2) 



/" ^ng ~^ uiuo"^ + iguiduo] = J dx 



pv 



+ P<P) 



(B.3) 



Higher integrals of motion for 2B0 can be constructed recurrently similarly to Benjamin- 
Ono equation [2S] or can be written using the integrals obtained for INLS 



Appendix C. Geometrical interpretation of the chiral equation: Contour 
dynamics 



Equation (67) can be cast in the form of contour dynamics. 

Let us interpret the unit disk as a uniformization of a simply-connected domain 
embedded into the complex z-plane. In other words, z{w) is a conformal map of the 
interior of the unit disk |w| = 1 to a bounded domain such that the length element of 
its boundary is proportional to the density 

ds = \z'{x)\dx = np{x)dx. (C.l) 



(equivalently 1/y^ is a harmonic measure of the contour). Then the equation (67) 
describes the evolution of the planar domain. We notice that the curvature of the 
boundary k = —ids log Zg can be expressed in terms of the density p{x) as 

K=-{7TprUdx{\0g^f-l). (C.2) 



Then (67) can be written as 



where d/dt = dt — gupod^ and the time derivative is taken at fixed x. The equation 



(C.3) describes the evolution of a planar contour driven by its curvature. 



CONTENTS 

Appendix D. Determinant identity 



25 



Here we present a derivation of a determinant identity (D.7) wliicli was used in Sec. 7.2 
Consider tlie Caucliy matrix 



Pi - Qj ^ 



PiT^Qj, = 1,2, ...,iV. 



Its determinant 



det(D) = Yl 



Ui<JiP^-Pj)iQ3-Q^ 



n'.= n 



fiSiTiSi I , 



where rj and Sj are defined by 



TiTi = {pi - qi) TT —, 

Pi — Pk 

kik^iy"- 



- Pi) n 

k(k^i) 



Qi-Pk 
Qi-Qk 



The inverse of D is also Cauchy matrix (see e.g., pS]) given by 



^j^i 

Pj - Qi 



The obvious identity 

det(l + D)= det{D) det(l + {D'^)'^] 
being specialized for Cauchy matrix D ( D.l[ ) reads 



iUi riSi) 



1/2 



^ S ^1/2 ■ 



(D.l) 
(D.2) 

(D.3) 
(D.4) 

(D.5) 
(D.6) 

(D.7) 



Appendix E. Non-degeneracy condition of matrix (107) 



In this section we will derive the non-degeneracy condition for the matrix (107) following 



Let us introduce 



Cj = cij exp 



9 



Our aim is to derive a sufficient condition on coefficients cj for the matrix 



Pj- Qk 

to be non-degenerate for real pj and qj. 



(E.l) 



(E.2) 



Let us assume that the matrix (E.2) is degenerate. It means the the following 
equation has a non-zero solution 



+ V — '—Tk = 0. 



(E.3) 
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We introduce the meromorphic function 



z-Qk 



and rewrite the condition (E.3) as 
Now, let us consider the function 



f{z) = ij{z)mll — 

Z — Pj 
3 ■' 



Here explicitly 



m = J2 



Tk 



Z-Qk 



The function f{z) is meromorphic and its residue at infinity is zero (as ip{z) 
z — >■ 06). On the other hand {pj and qj are real numbers) 

Pj - Qk 



1 T-r Qj-Qk _ b 



k{k^j) 



n 



^^-p^ktk;,)^^-^'^ 



.12 •'J' 



where we introduced a notation 



(E.4) 

(E.5) 
(E.6) 

(E.7) 

z""*" as 



(E.8) 



fj = (pj-Qj) n 



k(k^j) 



Pj - Pk 



1 - 



{Pj - Pk){qj - Qk) 
(Pj - (Ij? - (lk){(lj -Pk) 



n 



.(E.9) 



If fj has the same sign for all j, e.g., fj < we immediately obtain 



E (res^=p^/ + < 0. 



This is impossible as the residue of f{z) at infinity is zero. This contradiction shows 
that the matrix (E.2) is nondegenerate under these conditions. 



To summarize, the sufficient condition of non-degeneracy of (E.2) is 
fj < for all i or 



/, > for all J, 



(E.IO) 



where fj are defined according to (E.9). 
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